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The symmetry energy and the isoscaling properties of the fragments produced in the multifrag- 
mentation of 40 Ar, 40 Ca + 58 Fe, 58 Ni reactions at 25 - 53 MeV/nucleon were investigated within the 
framework of statistical multifragmentation model. The isoscaling parameters a, from the primary 
(hot) and secondary (cold) fragment yield distributions, were studied as a function of excitation 
energy, isospin (neutron-to-proton asymmetry) and fragment symmetry energy. It is observed that 
the isoscaling parameter a decreases with increasing excitation energy and decreasing symmetry 
energy. The parameter a is also observed to increase with increasing difference in the isospin of 
the fragmenting system. The sequential decay of the primary fragments into secondary fragments, 
when studied as a function of excitation energy and isospin of the fragmenting system, show very 
little influence on the isoscaling parameter. The symmetry energy however, has a strong influence 
on the isospin properties of the hot fragments. The experimentally observed scaling parameters can 
be explained by symmetry energy that is significantly lower than that for the ground state nuclei 
near saturation density. The results indicate that the properties of hot nuclei at excitation energies, 
densities and isospin away from the normal ground state nuclei could be significantly different. 
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I. INTRODUCTION 

Recently the possibility of extracting information on 
the symmetry energy and the isospin (neutron-to-proton 
ratio) of the fragments in a multifragmentation reaction 
has gained tremendous importance Q, 0] Such 
information is of importance for understanding some of 
the key problems in astrophysics 0, 0, II H EH EJ EI , 
and various aspects of nuclear physics such as the struc- 
ture of exotic nuclei (the binding energy and rms radii) 
US, EjL EM El and the dynamics of heavy ion collisions 
UA SMUfl IH El IM Ei- Traditionally, the sym- 
metry energy of nuclei has been extracted by fitting the 
binding energy in their ground state with various ver- 
sions of the liquid drop mass formula |28| . The proper- 
ties of nuclear matter are then determined by theoreti- 
cally extrapolating the nuclear models designed to study 
the structure of real nuclei. However, real nuclei are cold, 
nearly symmetric (N w Z) and found at equilibrium den- 
sity. It is not known how the symmetry energy behaves at 
temperatures, isospin (neutron-to-proton ratio) and den- 
sities away from the normal nuclear matter. Theoretical 
many-body calculations p9l l3fl l3ll l32| and those from 
the empirical liquid drop mass formula |33l l34| predict 
symmetry energy near normal nuclear density (« 0.17 
fm~ 3 ) and temperature (T w MeV), to be around 28 - 
32 MeV. 

In a multifragmentation reaction, an excited nucleus 
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expands to a sub-nuclear density and disi nteg rates into 
various light and heavy fragments [33, |3g, |33, |33| ■ The 
fragments are highly excited and neutron-rich; their 
yields depend on the available free energy which in turn 
depends on the strength of the symmetry energy and the 
extent to which the fragments expand. By studying the 
isotopic yield distribution of these fragments, one can ex- 
tract important information about the symmetry energy 
and the properties of the fragments at densities, excita- 
tion energies and isospin away from those of ground state 
nuclei. 

Experimentally, the determination of fragment isotopic 
yield distribution is not straight forward. It is influenced 
by the complex de-excitation of the hot (primary) frag- 
ments into observed cold (secondary) fragments. The- 
oretical calculations of the secondary de-excitations re- 
quire accurate accounting of the feeding from the parti- 
cle unstable states and are subject to uncertainties in the 
levels that can be excited, and the structure effects that 
govern their decay [3!J |4(|, ElJ . 

From various statistical model approaches 0, E^|, it 
has been shown that the ratio of primary fragment yield 
for a given isotope or isotone produced in two different 
reactions with similar temperature, exhibit an exponen- 
tial dependence on proton and neutron number, an ob- 
servation known as isoscaling. The dependence has been 
interpreted in terms of a scaling parameter that is related 
to the symmetry energy of the primary fragment bind- 
ing energy. The scaling parameter has been shown to be 
independent of the complex nature of the secondary de- 
excitation and is thus a robust observable for studying 
the fragment isotopic yield distribution. 
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In a recent work [l|, it was shown that the symme- 
try energy of the primary fragments deduced from the 
reduced neutron density is significantly lower than that 
for the normal nuclei at saturation density. Le Fevre 
et al. 2], in their recent work studied the fragmen- 
tation of excited target residues following collisions of 
12 C on 112 < 124 Sn, at incident energies of 300 and 600 
MeV/nucleon. They observed that the symmetry energy 
co-efficient deduced from the data are near 25 MeV for 
peripheral collisions and lower than 15 MeV for the cen- 
tral collisions. Hcnzlova et al. Q, studied the fragments 
produced in the multifragmentation of 136 ' 124 Xe projec- 
tiles in mid-peripheral collisions with a lead target at 1 
GeV/nucleon. They used both, the < N > jZ ratio and 
the isoscaling of the fragments and found that the exper- 
imentally determined value of the scaling parameter can 
be reproduced within the statistical model framework by 
lowering the symmetry energy to as low as 11 - 12 MeV. 
The < N > jZ ratios of the fragments on the other hand, 
could be reproduced with symmetry energy co-efficient as 
low as 14 - 15 MeV. 

In this work, we study the primary fragment yield 
distribution in a number of reactions at various ex- 
citation energies and isospins using the isoscaling ap- 
proach and the equilibrium statistical multifragmenta- 
tion model. It is observed that the isoscaling parameter a 
for the hot fragments decreases with increasing excitation 
energy and decreasing symmetry energy. The alpha val- 
ues also increase with increasing difference in the isospin 
(neutron-to-proton asymmetry) of the fragmenting sys- 
tem. A similar behavior is also observed for the cold sec- 
ondary fragments. The secondary de-excitation is found 
to have very little influence on the isoscaling parameter at 
lower excitation energies and isospin of the fragmenting 
system. The symmetry energy however, strongly influ- 
ence the properties of the hot fragments. The experimen- 
tally determined isoscaling parameters can be explained 
by symmetry energy that is significantly lower than that 
for the normal (cold) nuclei at saturation density, indi- 
cating that the properties of nuclei at high excitation 
energy, isospin and reduced density are sensitive to the 
symmetry energy. 

The paper is organized as follows. In sec. II, we de- 
scribe the experiment in detail. The experimental results 
are presented and discussed in Sec. III. Section IV con- 
tains a brief description of the statistical multifragmen- 
tation model used in the analysis. Section V contains a 
comparison between the experimental data and the sta- 
tistical multifragmentation model. Finally, a summary 
and conclusion is given in section VI. 



II. EXPERIMENT 

A. Experimental Setup 

The experiments were carried out at the Cyclotron 
Institute of Texas A&M University (TAMU) using the 




FIG. 1: Schematic diagram of the experimental setup showing 
the placement of the telescopes inside the scattering chamber 
of the neutron ball. 



K500 Superconducting Cyclotron where targets of 58 Fe 
(2.3 mg/cm 2 ) and 58 Ni (1.75 mg/cm 2 ) were bombarded 
with beams of 40 Ar and 40 Ca at 25, 33, 45 and 53 
MeV/nucleon [4^|. The various combinations of target 
and projectile nuclei allowed for a range of N/Z (neutron- 
to-proton ratio) (1.04 - 1.23) of the system to be stud- 
ied, while keeping the total mass constant (A = 98). 
The beams were fully stripped by allowing them to pass 
through a thin aluminum foil before being hit at the cen- 
ter of the target inside the TAMU 47r neutron ball [43| . 
Light charged particles (Z < 2) and intermediate mass 
fragments (Z > 2) were detected using six discrete tele- 
scopes placed inside the scattering chamber of the neu- 
tron ball at angles of 10°, 44°, 72°, 100°, 128° and 148°. 
A schematic diagram of the setup showing placement of 
the telescopes within the scattering chamber of the neu- 
tron ball is shown in figure. 1. Each telescope consisted 
of a gas ionization chamber (IC) followed by a pair of 
silicon detectors (Si-Si) and a Csl scintillator detector, 
providing three distinct detector pairs (IC-Si, Si-Si, and 
Si-Csl) for fragment identification. The ionization cham- 
ber was of axial field design and was operated with CF4 
gas at a pressure of 50 Torr. The gaseous medium was 6 
cm thick with a typical threshold of ~ 0.5 MeV/nucleon 
for intermediate mass fragments. The silicon detectors 
had an active area of 5 cm x 5 cm and were each sub- 
divided into four quadrants. The first and second silicon 
detectors in the stack were 0.14 mm and 1 mm thick, 
respectively. The dynamical energy range of the silicon 
pair was ~ 16 - 50 MeV for 4 He and - 90 - 270 MeV 
for 12 C. The Csl scintillator crystals that followed the 
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silicon detector pair were 2.54 cm in thickness and were 
read out by photodiodes. Good elemental (Z) identifi- 
cation was achieved for fragments that punched through 
the IC detector and stopped in the first silicon detector. 
Fragments measured in the Si-Si detector pair also had 
good isotopic separation. The trigger for the data acqui- 
sition was generated by requiring a valid hit in one of the 
silicon detectors. 



B. Calibration 

The calibration of the IC-Si detectors were carried out 
using the standard alpha sources and by operating the 
IC at various gas pressures. The Si-Si detectors were cal- 
ibrated by measuring the energy deposited by the alpha 
particles in the thin silicon and the punch-through ener- 
gies of different isotopes in the thick silicon. The Si-Csl 
detectors were calibrated by selecting points along the 
different light charged isotopes and determining the en- 
ergy deposited in the CSI crystal from the energy loss in 
the calibrated Si detector. 



C. Neutron Multiplicity and Event 
Characterization 

Neutrons were measured with the 47r neutron ball that 
surrounded the detector assembly. The neutron ball con- 
sisted of eleven scintillator tanks segmented in its median 
plane and surrounding the vacuum chamber. The upper 
and the lower tank were 1.5 m diameter hemispheres. 
Nine wedge-shaped detectors were sandwiched between 
the hemispheres. All the wedges subtended 40° in the 
horizontal plane. The neutron ball were filled with a 
pscudocumene based liquid scintillator mixed with 0.3 
weight percent of Gd salt (Gd 2-ethyl hexanoate). Scin- 
tillations from thermal neutrons captured by Gd were 
detected by twenty 5-in phototubes : five in each hemi- 
sphere, one on each of the identical 40° wedges and two 
on the forward edges. The efficiency with which the neu- 
trons could be detected is about 83%, as measured with 
a 252 Cf source. 

The detected neutrons were used to differentiate be- 
tween the central and peripheral collisions. To un- 
derstand the effectiveness of neutron multiplicity as 
a centrality trigger, simulations were carried out us- 
ing a hybrid BUU-GEMINI calculations at various im- 
pact parameters for the 40 Ca + 58 Fe reaction at 33 
MeV/nucleon. The simulated neutron multiplicity distri- 
bution was compared with the experimentally measured 
distribution. The multiplicity of neutron for the impact 
parameter 6 = collisions was found to be higher than 
the 6 = 5 collision. By gating on the 10% highest neu- 
tron multiplicity events, one could clearly discriminate 
against the peripheral events. 

To determine the contributions from noncentral im- 
pact parameter collisions, neutrons emitted in coinci- 
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FIG. 2: Relative isotopic yield distribution for the Lithium 
(left), Berillium(center) and Carbon elements in 40 Ca + 58 Ni 
(stars and solid lines), 40 Ar + 58 Ni (circles and dashed lines) 
and 40 Ar + 58 Fe (squares and dotted lines) reactions at vari- 
ous beam energies. 



dence with fragments at 44° and 152° were calculated at 
6 = fm and 6 = 5 fm. The number of events were ad- 
justed for geometrical cross sectional differences. A ratio 
was made between the number of events with a neutron 
multiplicity of at least six, calculated at b = fm, and the 
number of events with the same neutron multiplicity at b 
= 5 fm. The ratios were observed to be 19.0 and 11.1 at 
44° and 1.3 and 2.2 at 152° for 33 and 45 MeV/nucleon 
respectively. At intermediate angles, high neutron multi- 
plicities were observed to be outside the region in which 
6 = 5 fm contributes significantly. At backward angles 
the collisions at 6 = 5 fm made a larger contribution to 
the neutron multiplicity. 

In addition to the neutron multiplicity distribution, the 
charge distribution of the fragments was also used to in- 
vestigate the contributions from central and mid-impact 
parameter collisions. The 6 = 5 collisions produced es- 
sentially no fragments with charge greater than three in 
the 44° telescope. 

In an earlier work 42§], some analysis of the fragment 
kinetic energy and charge distributions were presented. 



4 



It was shown that at a laboratory angle of 44° the ki- 
netic energy and the charge distributions are well repro- 
duced by the statistical model calculation. Using a mov- 
ing source analysis of the fragment energy spectra, it was 
also shown that the fragments emitted at backward an- 
gles originate from a target-like source, while those emit- 
ted at 44° originate primarily from a composite source. 
In this work, we will concentrate exclusively on data from 
the laboratory angle of 44° to study the symmetry energy 
and the isoscaling properties of the fragments produced. 
The choice of this angle enables one to select events which 
are predominantly central and undergo bulk multifrag- 
mentation. The contributions to the intermediate mass 
fragments from the projectile-like and target-like sources 
can therefore be assumed to be minimal and the use of 
equilibrium statistical model appropriately justified. 



III. EXPERIMENTAL RESULTS 

A. Relative Fragment Yield 

The experimentally measured relative isotopic yield 
distributions for the Lithium (left), Berillium (center) 
and Carbon (right) elements, in 40 Ca + 58 Ni (star sym- 
bols), 40 Ar + 58 Ni (circle symbols) and 40 Ar + 58 Fe 
(square symbols) reactions, are shown in figure 2 for the 
beam energies of 25, 33 and 45 MeV/nucleon. The dis- 
tributions for each element show higher fragment yield 
for the neutron rich isotopes in 40 Ar + 58 Fe reaction 
(squares), which has the largest neutron-to-proton ra- 
tio (N/Z), in comparison to the 40 Ca + 58 Ni reaction 
(stars), which has the smallest neutron-to-proton ratio. 
The yields for the reaction, 40 Ar + 58 Fe (circles), which 
has an intermediate value of the neutron-to-proton ra- 
tio, are in between those of the other two reactions. The 
figure thus shows the isospin dependence of the compos- 
ite system on the properties of the fragments produced 
in the multifragmentation reaction. One also observes 
that the relative difference in the yield distribution be- 
tween the three reactions decreases with increasing beam 
energy. This is due to the secondary de-excitation of 
the primary fragments, a process that becomes impor- 
tant for systems with increasing neutron-to-proton ratio 
and excitation energy. In the following subsections, we 
will utilize the experimentally determined isotopic yield 
distributions to establish the isoscaling properties of the 
produced fragments before comparing them with the sta- 
tistical multifragmentation model in section V. 



B. Isotopic and Isotonic Scaling 

In a multifragmentation reactions, the ratio of isotope 
yields in two different systems, 1 and 2, R 21 (N 7 Z) = 
Y 2 (N ', Z) I 'Y 1 (N ', Z) , has been shown to obey an expo- 
nential dependence on the neutron number (N) and the 
proton number (Z) of the isotopes, an observation known 
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FIG. 3: Double isotope ratio as a function of the difference in 
binding energy for various beam energies. The circle symbols 
correspond to 40 Ca + 58 Ni reaction, square symbols to 40 Ar + 
58 Ni reaction, and triangle symbols to 40 Ar + 58 Fe reaction. 
The solid lines are the best fit to the data. 



as isoscaling 0, |2|| 0] . The dependence is charac- 
terized by a simple relation, 



R 21 (N,Z)=Y 2 (N,Z)/Y 1 (N,Z) = Cexp{aN + 0Z) 

(1) 

where, Y 2 and Y\ are the yields from the neutron- 
rich and neutron-deficient systems, respectively. C, is 
an overall normalization factor and a and are the pa- 
rameters characterizing the isoscaling behavior. 

The necessary condition for observing isoscaling in 
multifragmentation reaction is the near equality of tem- 
perature for systems chosen in the above scaling re- 
lation [451 ]. For the present work, this condition was 
tested by determining the temperature usin g th e dou- 
ble isotope ratio method of Albergo et al. 46] . The 
Albergo method relates the apparent temperature T of 
the system at the freeze-out to the double isotope ratio, 
R = [Y(N, Z 1 )/Y(N + 1, Z 1 )]/[Y(N, Z 2 )/Y{N + 1, Z 2 )], 
through a relation 
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FIG. 4: Experimental isotopic yield ratios of the fragments as 
a function of neutron number N, for various beam energies. 
The left column correspond to 40 Ar + 5S Fe and 40 Ca + 5S Ni 
pair of reactions. The right column correspond to 40 Ar + 58 Ni 
and 40 Ca + 58 Ni pair of reactions. The different symbols 
correspond to Z — 3 (circles), Z = 4 (open stars), Z = 5 
(triangles), Z = 6 (squares) and Z = 7 (filled stars) elements. 
The lines are the exponential fits to the data as explained in 
the text. 



FIG. 5: Experimental isotonic yield ratios of the fragments 
as a function of proton number Z, for various beam energies. 
The left column correspond to 40 Ar + 5S Fe and 40 Ca + 58 Ni 
pair of reaction. The right column correspond to 40 Ar + 58 Ni 
and 40 Ca + 58 Ni pair of reaction. The different symbols cor- 
respond to iV = 3 (circles) , N = 5 (triangles) , N = 6 (squares) 
and N = 7 (stars) elements. The lines are the exponential fits 
to the data as explained in the text. 
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where a is a factor that depends on the statistical 
weights of the ground state nuclear spins, and AB = 
[B(N, Zi) - B(N + 1, Zi)] - [B(N, Z 2 ) - B(N + 1, Z 2 )] 
is the difference in the binding energy. The method was 
applied to all the three systems and beam energies stud- 
ied. 

Figure 3 shows the double isotope ratios obtained us- 
ing the various combinations of isotopes as a function 
of the difference in the binding energy. The top panel 
correspond to those determined for the beam energies 
of 25 MeV/nucleon, the center for the 45 MeV/nucleon 
and the bottom for the 53 MeV/nucleon. The differ- 



ent symbols correspond to the three different reactions 
studied. The isotope ratios for the three reactions at 
each beam energy overlap nicely indicating formation of 
composite systems with similar temperature/excitation 
energy. The slopes of the exponential fits to the data 
(shown by solid line), which correspond to the apparent 
temperatures of the systems, show a gradual decrease 
with increasing beam energy indicating an increase in 
the excitation energy/temperature of the system. 

Having satisfied the necessary condition for isoscal- 
ing, the isotopic yield ratio as a function of neutron 
number N, for the beam energies of 25, 33, 45 and 53 
MeV/nucleon is plotted in Fig 4. The left column shows 
the ratio for the 40 Ar + 58 Fe and 40 Ca + 58 Ni pair of re- 
action and the right column shows the ratio for the 40 Ar 
+ 58 Ni and 40 Ca -I- 58 Ni pair of reaction. One observes 
that the ratio for each element lies along a straight line 
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FIG. 6: Isoscaling parameters a (solid symbols) and /3 (open 
symbols) as a function of the beam energy. The solid circles 
and open squares correspond to 40 Ar + 58 Fe and 40 Ca + 
58 Ni pair of reactions. The solid stars and open triangles 
correspond to 40 Ar + 58 Ni and 40 Ca + 58 Ni pair of reactions. 
The lines are the exponential fits to the data. The error bars 
are of the size of the symbols. 



in the logarithmic plot and align with the neighboring 
element quite well. This feature is observed for all the 
beam energies and both pairs of reaction studied. One 
observes that the alignment of the data points varies with 
beam energies as well as the pairs of reaction. To have a 
quantitative estimate of this variation, the ratio for each 
elements (Z) were simultaneously fitted using an expo- 
nential relation (shown by the solid lines) to obtain the 
slope parameter a. The values of the parameters are 
shown at the top of each panel in the figure. The value 
of the slope parameter a is larger for the 40 Ar + 58 Fe and 
40 Ca + 58 Ni reactions, which has a larger difference in 
the N/Z of the systems in the pair, compared to the 40 Ar 
+ 58 Ni and 40 Ca + 58 Ni reactions, which has a smaller 
difference in the corresponding N/Z. The a value fur- 
thermore decreases with increasing beam energy. Figure 
5. shows the isotonic yield ratio as a function of atomic 
number Z, for the same beam energies and pairs of sys- 
tems as shown in Figure 4. Once again, one observes the 
ratios for each isobar to align nicely with each other at 
all beam energies. The scaling parameters j3 in this case, 
shows an increase with increasing beam energy. The val- 
ues of the j3 parameter are larger for the 40 Ar + 58 Ni and 
40 Ca + 58 Ni reaction pair compared to the 40 Ar + 58 Fe 
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The top panel is for the 40 Ar + 58 Fe/ 40 Ca + 58 Ni 



pair, and the bottom panel is for the °Ar + 58 Ni/ 4U Ca + 
58 Ni pair of reactions. The symbols correspond to S(f5,N) 
obtained from various elements (Z). The lines are the best 
fits to the data. 



and 40 Ca + 58 Ni reaction pair. Figure 6 shows a relative 
comparison of how the a and j3 parameters evolve as a 
function of beam energy and the isospin of the system. 

The temperature/excitation energy dependence of the 
isoscaling properties can be further studied by construct- 
ing the scaling factor, S(f3,N) = R 21 (N, Z)e~P z . The 
scaling factor is known to be a robust feature over a 
large range of data, from deep inelastic heavy ion re- 
actions at lowest energies, through evaporation reactions 
induced by light-ion and heavy-ion projectiles, to high 
energy heavy ion reactions characterized by intermedi- 
ate mass fragments and multifragmentation j2(J . With a 
single value for the j3 parameter, all the isotopes should 
fall along a single line in a plot of S((3, N) vs N. This 
is shown in Fig. 7, where the S(/3, N) from two differ- 
ent reaction pairs are plotted as a function of neutron 
number N for beam energies of 25 MeV/nucleon and 45 
MeV/nucleon. The parameter /?, was taken from the fit 
to the isotonic yield ratio shown in figure 5. As shown 
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in figure 7, the values of S(f3, N) obtained from various 
elements (Z) cluster and scale along a single line. The 
figure shows a significant difference in the scaling for the 
two beam energies, indicating the influence of temper- 
ature on the isotopic yields of the light clusters. The 
present observation alternatively demonstrates the role 
played by the temperature in the distillation of nuclear 
matter into a neutron-rich gas and a symmetric liquid 
phase. We illustrate this further in the following section. 



Isospin fractionation and the reduced nucleon 
densities 



In Grand-Canonic al appro ach of the multifragmenta- 
tion process (see e.g. |i(||42],|43,|4{|), the fragment yield 
with neutron number N, and proton number Z (mass 
number A = N + Z), can be written as 



Y(N, Z) oc Vp^p^Z N ,z{T)A^e B ^l T (3) 

where V is the volume of the system and p n (oc e M ™/ T ) 
and p p (oc e^ p / T ) are the primary 'free' neutron and pro- 
ton densities. The exponents p n and p p are the neutron 
and the proton chemical potentials, and Z^^z{T) is the 
intrinsic partition function of the excited fragment. The 
quantity B(N, Z), is the ground state binding energy of 



the fragment and T is the temperature. In the above 
formula, the effect of Coulomb interaction on fragment 
yield is neglected by introducing p n and p p . The actual 
isotope yields then reduce to an approximation appro- 
priate for the thermodynamical limit at high excitation 
energy Jijj]. As discussed in the introduction, taking the 
ratios of the fragment yields from two different systems 
which differ only in their isospin (N/Z) content reduces 
the uncertainties in the quantities shown in Eq. 3. The 
isotopic yield distribution of the fragments in terms of 
relative reduced neutron density can then be written as, 



Y(N + k, Z) /Y Ca+m (N + k,Z) _ 
YjN, Z)/Y Ca + Nl {N, Z) ~ 



Pn 

Ca+Ni 
Pn 



(4) 



where k, corresponds to various isotopes of an element 
that can be used to determine the double ratio. The 
quantity Y Ca+Ni is the yield for the 40 Ca + 58 Ni re- 
action with respect to which all the ratios are taken in 
this work. A similar expression for the relative reduced 
proton density from the isotonic yield ratios can also be 
written as 



Y(N, Z + k)/Y Ni+m (N, Z + k) 
Y(N, Z)/Y Ca+m {N, Z) 



Ca+Ni 
Pp 



(5) 



In the statistical limit for a dilute non-interacting gas, 
the relative nucleon densities are related to the isoscaling 
parameters a and /3, through a relation, Pn T+Fe / Pn a+N% 
= e a and p£ r+Fe / p^ a+Nl = e?, where a = Ap n /T and 
/3 = App/T, with Ap n and Ap p being the difference in 
neutron and proton chemical potentials. 

In figure 8, we show the experimentally obtained rel- 
ative reduced neutron and proton density as a function 
of the difference in the N/Z of the systems for the Ar + 
Ni/Ca + Ni and Ar + Fc/Ca + Ni pairs of reactions. All 
densities shown are relative to those of Ca + Ni reaction. 
The circle symbols correspond to the 25 MeV/nucleon, 
squares to the 33 MeV/nucleon and the triangles to the 
45 MeV/nucleon beam energies. The a and the (3 values 
were taken from the fit to the isotopic and isotonic yield 
ratios of figures 4 and 5. The figure shows a steady de- 
crease in the reduced neutron density and an increase in 
the proton density with increasing beam energies. The 
effect is stronger for the 40 Ar + 58 Fe and 40 Ca + 58 Ni 
reaction pair which has the highest difference in N/Z. 

An important feature of the data shown in figure 8, 
is the decrease in the relative neutron-proton asymmetry 
(p n / Pn a+Nl ~ Pp/ Pp a+Nl ) with increasing beam energy. 
The asymmetry is found to decrease from ~ 1.0 at 25 
MeV/nucleon to — 0.6 at 45 MeV/nucleon for the Ar + 
Fe/Ca -I- Ni pair of reaction. This is consistent with the 
prediction of phase coexistence of two conserved charges 
with isospin dependent interaction of lattice gas calcu- 
lation [53, |54j . The large asymmetry or the strong en- 
richment of the gas phase can be understood in terms 
of an isospin dependent interaction between the protons 
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and the neutrons. In the absence of an isospin depen- 
dent interaction, the fragments are populated close to 
the stability line and have similar N/Z ratio's as the ini- 
tial source. With the inclusion of isospin dependence, 
the heavy fragments are strongly favored from the ener- 
getic point of view to be populated closer to the bottom 
of the stability valley. Only a few protons are left to 
be shared between the light fragments leading to the ob- 
served enhancement of neutrons in the gas phase. With 
the increase in temperature, the fragment distribution is 
centered closer to the N = Z line resulting in a decrease 
in the asymmetry of the gas phase. 

The observed decrease in the relative neutron and pro- 
ton densities with increasing beam energy may thus be 
attributed to the decrease in the sensitivity of the isospin 
effect with increasing temperature. The scaling parame- 
ter is also sensitive to the breakup density, and the ob- 
served decrease in the neutron-proton density could re- 
sult from a decrease in the break-up density, though it is 
expected to be less sensitive. 



E A,z = l( A - 2Z ) 2 / A i where 7 = 25 MeV is the sym- 
metry energy co-efficient. These parameters are those 
adopted from the Bethe-Weizsacker mass formula and 
correspond to the assumption of isolated fragments with 
normal density in the freeze-out configuration. The value 
of the symmetry energy co-efficient 7 is taken from the 
fit to the binding energies of isolated cold nuclei in their 
ground states. In a multifragmentation process the pri- 
mary fragments are not only excited but also expanded. 
The fragments continue to interact in-medium with each 
other inside the freeze-out volume and modify their pa- 
rameters. By comparing the experimentally determined 
fragment yield distribution with the SMM calculation, 
the parameters of hot nuclei under multifragmentation 
conditions, including the symmetry energy, can be ex- 
tracted. In the following, it will be shown how this infor- 
mation can be obtained from the isoscaling phenomena. 



A. Isoscaling and Symmetry Energy Coefficient 



IV. STATISTICAL MULTIFRAGMENTATION 
MODEL 

Statistical models HI El |U E3, EJ E3 are widely 
used for describing multifragmentation reactions [321 IH, 
M, E3, ES El El They are based on the assump- 
tion of statistical equilibrium at a low density freeze- 
out stage. In the Statistical Multifragmentation Model 
(SMM) 0, E3 1 all break-up channels composed of nu- 
cleons and excited fragments are taken into account and 
considered as partitions. During each partition, the con- 
servation of mass, charge, energy and angular momen- 
tum is taken into account, and the partitions are sam- 
pled uniformly in the phase space according to their sta- 
tistical weights using the Monte Carlo sampling. In the 
present calculations the Coulomb interaction between the 
fragments is treated in the Wigner-Seitz approximation. 
Light fragments with mass number A < 4 are considered 
as elementary particles with only translational degrees 
of freedom ("nuclear gas"). Fragments with A > 4 are 
treated as heated nuclear liquid drops, and their individ- 
ual free energies Fa,z are parametrized as a sum of the 
volume, surface, Coulomb, and symmetry energy, 

Fa,z = FX >Z + F s AtZ + E% z + E%2 (6) 

where F\ z = (—W — T 2 /e )A 7 with parameter e a 
related to the level density and W = 16 MeV being 
the binding energy of infinite nuclear matter. F^ z — 

B A 2 / 3 [(T 2 - T 2 )/(T 2 + T 2 )] 5 / 4 , with B = 18 MeV 
being the surface co-efficient and T c = 18 MeV be- 
ing the critical temperature of infinite nuclear matter. 
E% z = cZ 2 /AV 3 , where c= (3/5)ie?/r )[l-(p/Po) 1/s \, 
is the Coulomb parameter obtained in the Wigner-Seitz 
approximation with charge unit e, and r Q — 1.17 fm. 



Isotopic scaling or isoscaling arise naturally in statis- 
tical equilibrium models of multifragmentation. In these 
models the difference in the chemical potential of sys- 
tems with different N/Z is directly related to the scaling 
parameter a. It has been shown that the isoscaling pa- 
rameter a is proportional to the symmetry energy part 
of the fragment binding energy through a relation, 



47 (zl zi\ 

where Z±, A\ and Z2, A2 are the charge and the mass 
numbers of the fragmenting systems, T is the temper- 
ature of the system and 7, the symmetry energy co- 
efficient 0- 



B. Secondary de-excitation of the fragments 

The above formula in the statistical model approach 
is valid at the freeze-out stage where the primary hot 
fragments are formed at reduced density. In order to 
extract information on symmetry energy 7, from the ob- 
served cold secondary fragments, one has to take into ac- 
count the process of secondary de-excitation. In SMM, 
the secondary de-excitation of large fragments with A 
> 16 is described by Weisskopf-type evaporation and 
Bohr- Wheeler- type fission models (471 |49j. The decay 
of smaller fragments is treated with the Fermi-breakup 
model. All ground and nucleon-stable excited states of 
light fragments are taken into account and the popula- 
tion probabilities of these states are calculated according 
to the available phase space [49j . 
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FIG. 9: SMM calculated primary (left) and secondary (right) 
fragment isotope yield distributions for the carbon element 
in 40 Ca + 58 Ni (stars and solid lines), 40 Ar + 58 Ni (circles 
and dashed lines) and 40 Ar + 58 Fe (squares and dotted lines) 
reactions at various excitation energies. The calculations are 
for 7 = 25 MeV. 



V. COMPARISON WITH THE 
EXPERIMENTAL DATA 



In order to compare the experimentally observed re- 
sults to the theoretical predictions, we have calculated 
the primary and the secondary fragment isotopic yield 
distributions using the above described SMM. The calcu- 
lations were carried out for the 40 Ca + 58 Ni, 40 Ar + 58 Ni 
and 40 Ar + 58 Fe reactions at various excitation energies. 
The excitation energy per nucleon of the initial system 
for the calculation depends strongly on the matching con- 
dition between the dynamical and statistical stage of the 
collision. This quantity is presently difficult to calculate 
accurately. A range of values for the excitation energy 
per nucleon from E* = 4 - 10 MeV/nucleon were there- 
fore assumed. The excitation energy corresponding to 
each beam energy was also verified by an independent cal- 
culation using BUU - GEMINI (see table II in Ref.|l|), 
and the systematic calorimetric measurements available 
in the literature for systems with mass A ~ 100 |6l|. The 



FIG. 10: Calculated isotopic yield ratios from the primary 
(left) and the secondary (right) fragment yield distributions 
for the 40 Ar + 58 Ni and 40 Ca + 58 Ni pairs at various exci- 
tation energies. The calculations are for 7 = 25 MeV. The 
different symbols shown correspond to Z — 3 (circles), Z = 
4 (open stars), Z = 5 (triangles), Z — 6 (square) and Z = 
7 (solid stars) elements. The lines are the exponential fits to 
the ratios. 



mass and the charge of the initial systems were assumed 
to be those of the initial compound nucleus. To check 
for the possible uncertainty in the source size due to loss 
of nucleons during pre-equilibrium emission, the calcula- 
tion was also performed for sources with 80% of the total 
mass. No significant change in the isospin characteristics 
under study was observed. The freeze-out density in the 
calculation was assumed to be 1/3 of the normal nuclear 
density and the symmetry energy co-efficient 7 was taken 
to be 25 MeV. 

The calculated primary and secondary fragment yield 
distributions for the Carbon isotopes in 40 Ca + 58 Ni, 
40 Ar + 58 Ni and 40 Ar + 58 Fe reactions at various ex- 
citation energies are as shown in figure 9. The charac- 
teristics of the hot primary fragment yield distribution, 
shown in the left column of the figure, change signifi- 
cantly after the secondary de-excitation. The primary 
yield distribution for the three systems, shown by the 
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FIG. 11: Comparison of the SMM calculated a (lines) with 
the experimentally determined a (symbols) as a function of 
excitation energy for different values of the symmetry energy 
co-efficient 7. The dotted lines correspond to the primary 
fragments and the solid lines to the secondary fragments. The 
left column shows the comparison for the 40 Ar + 58 Ni and 
40 Ca + 58 Ni pair, and the right column shows the comparison 
for the 40 Ar + 58 Fe and 40 Ca + 58 Ni pair. 



dotted, dashed and solid curve, for each excitation energy 
clearly show isospin effect. The most neutron rich sys- 
tem has the largest yield for the neutron rich isotopes and 
the lowest yield for the neutron deficient isotopes. This 
effect appears to decrease with increasing excitation en- 
ergy. A similar feature is also observed in the secondary 
fragment distribution shown in the right column of the 
figure, though the effect is observed to be weakened sig- 
nificantly. Furthermore, the mean of the distribution is 
also observed to decrease along with the width for the 
secondary fragments. Qualitatively, the SMM simulates 
quite well the overall features of the experimentally ob- 
served isotopic yield distribution shown in figure 2. The 
isotopic yield ratios using the primary and the secondary 
fragment distribution from the statistical multifragmen- 
tation model are shown in figure 10 for the 40 Ar + 58 Ni 
and 40 Ca + 58 Ni pair of reaction. It is observed that the 
yield ratios for both the primary and the secondary distri- 
bution obey the isoscaling relation quite well. Very little 
difference in the scaling parameter a, obtained from the 
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FIG. 12: SMM calculated isoscaling parameter a as a function 
of symmetry energy co-efficient for various excitation energies. 
The open circles joined by dotted lines correspond to the pri- 
mary fragments and the open stars joined by solid lines to the 
secondary fragments. The left column shows the calculation 
for 40 Ar + 58 Ni and 40 Ca + 58 Ni pair, and the right column 
for the 40 Ar + 58 Fe and 40 Ca + 58 Ni pair. 



primary and the secondary yield distribution is observed. 
Furthermore, the scaling parameter shows a gradual de- 
crease in its value with increasing excitation energy simi- 
lar to those observed experimentally and shown in figure 
6. 

Though the overall feature of the scaling parameter 
calculated from the statistical multifragmentation model 
is reproduced quite well, the absolute values do not quite 
agree with the experimentally determined a. This is 
shown in the top panel of figure 11, where a compari- 
son is made between the SMM calculated and the ex- 
perimentally observed values of a. In the figure, the left 
column correspond to the 40 Ar + 58 Ni and 40 Ca + 58 Ni 
pair of reaction and the right to the 40 Ar + 58 Fe and 
40 Ca + 58 Ni pair. The dotted lines correspond to a cal- 
culated from the primary fragment distribution and the 
solid lines to those calculated from the secondary frag- 
ment distribution. The symbols correspond to the ex- 
perimentally determined a's. It is observed that the ex- 
perimentally determined a's are significantly lower than 
the calculated values of a using the standard value of 
the symmetry energy co-efficient, 7 — 25 MeV, for the 
isolated cold nuclei in their ground states. In order to 
explain the observed dependence of the isoscaling param- 
eter a on excitation energy, we varied the 7 of the hot 
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FIG. 13: Same as in fig. 11, but with the modified secondary 
de-excitation with evolving symmetry energy co-efficient. 



primary fragment in the SMM input in the range 25 - 
15 MeV. As shown in the center and the bottom panel 
of the figure, the isoscaling parameter decreases slowly 
with decreasing symmetry energy. The experimentally 
determined a could be reproduced for both pairs of sys- 
tems at all excitation energies using a symmetry energy 
value of 7 = 15 MeV. This value of symmetry energy is 
significantly lower than the value of 7 = 25 MeV used for 
ground state nuclei. 

Figure 12. shows the calculated 7 dependence of the 
isoscaling parameter a from the hot primary fragment 
distribution, and from the cold secondary fragment dis- 
tribution, for the two pairs of systems at various excita- 
tion energies. The a as a function of symmetry energy 
for each excitation energy and system is observed to de- 
crease with decreasing symmetry energy. The difference 
between the primary fragment a and the secondary frag- 
ment a is negligible for the Ar + Ni and Ca + Ni reaction 
pair, which has the lowest difference in neutron-to-proton 
ratio. The difference for the Ar + Fe and Ca + Ni pair, 
which has the highest difference in neutron-to-proton ra- 
tio, however is slightly larger at higher excitation ener- 
gies. 



A. Secondary de-excitation of the fragments with 
changing symmetry 

In the above described statistical multifragmentation 
model calculations the masses of the fragments used were 
those of cold isolated nuclei. The fragments in their pri- 
mary stage are usually hot and the properties of hot nu- 
clei (i.e., their binding energies and masses) differ from 
those of cold nuclei. If hot fragments in the freeze-out 
configuration have smaller 7, their masses at the begin- 
ning of the secondary de-excitation will be different, and 
this effect should be taken into account in the evapora- 
tion process. Recently, Buyukcizmeci et al. |63| adopted 
a phcnomcnological approach to estimate the effect of the 
symmetry energy evolution during the sequential evapo- 
ration. In this approach, they assume liquid drop masses 
= mud) f° r the evaporation of the light particles 
(n, p, d, t, 3 He, a), if the internal excitation energy of 
the fragment is large (£ = f3E* /A > 1). At lower exci- 
tation energies (£ < 1) they assume a smooth transition 
to standard experimental masses with shell effects(m e2: p) 
using the following dependence, 



m A ,z = mid{l)£ + m exp (l - f) 



(8) 



The excitation energy is determined from the energy 
balance taking into account the mass uia,z at the given 
excitation. The above corrections were incorporated in 
the statistical model calculations described in the previ- 
ous section to study the effect of 7 during the sequential 
de-excitation of the hot primary fragments. Figure 13 
shows the result of the statistical model calculation us- 
ing the modified secondary de-excitation correction. The 
isoscaling parameter a is plotted as a function of excita- 
tion energy for the two pairs of systems (left and right 
column). The top panels show the calculations using 
symmetry energy value of 25 MeV. As noted in the pre- 
vious section, the new calculations are not able to repro- 
duce the experimentally determined alpha for both pairs 
of systems. With decreasing values of the symmetry en- 
ergy, the calculated a values for the Ar + Ni and Ca + Ni 
pair (left column) decrease and are in better agreement 
with the experimental values at 7 = 13 MeV. On the 
other hand, the calculated values for the Ar + Fe and Ca 
+ Ni reaction pair (right column) are in good agreement 
with the experimental values at 7 = 10 MeV. In gen- 
eral, one observes that the modified version of secondary 
de-excitation in SMM leads to a symmetry energy value 
of 10 - 13 MeV. This is slightly lower than the value of 
15 MeV obtained from the standard version of the SMM 
calculation shown in figure 11. 

The dependence of the isoscaling parameter as a func- 
tion of the symmetry energy for primary and secondary 
fragments at various excitation energies are as shown in 
figure 14. Once again, the difference between the primary 
fragment a and the secondary fragment a are extremely 
small for the Ar + Ni and Ca + Ni pair of reaction (which 
has the lowest difference in neutron-to-proton ratio), and 
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FIG. 14: Same as in fig. 12, but with the modified secondary 
de-excitation with evolving symmetry energy co-efficient. 



only slightly larger for the Ar + Fe and Ca + Ni pair 
(which has the highest difference in neutron-to-proton 
ratio). However, the main difference between the depen- 
dence shown in figure 12 and 14 is the rate at which the 
isoscaling parameter a decrease with decreasing symme- 
try energy. The decrease is much slower in the calcula- 
tion where the symmetry energy dependence of the mass 
is taken into account during the secondary de-excitation. 
The slower decrease in the isoscaling parameter results in 
the calculation being able to reproduce the experimental 
value at a slightly lower value of symmetry energy. One 
also observes that the primary fragment a in the modi- 
fied version of the secondary decay calculation are con- 
sistently smaller than the secondary fragment a at lower 
symmetry energies. 

In order to understand the difference in the symme- 
try energy, we show in figure 15, the calculated isotopic 
yield distribution for the Carbon element in 40 Ar + 58 Fe 
reaction at E* = 6 MeV/nucleon. The figure shows the 
primary and the secondary yield distribution using two 
different prescriptions for the secondary de-excitations in 
SMM. The left column in the figure corresponds to the 
SMM calculations, where the fragment masses used are 
those of cold isolated nuclei, and the right column corre- 
sponds to the SMM calculations, where the masses (sym- 
metry energy) evolve with their excitation energy during 
secondary de-excitation. The panels on the top corre- 
spond to the primary yield distribution and those in the 
bottom to the secondary yield distribution. The dotted 
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and the solid curves in each panel corresponds to the cal- 
culations assuming two different values of the symmetry 
energy, 15 and 25 MeV respectively. From the figure, it 
is evident that there exist a subtle difference between the 
two final (secondary) distributions, shown by the dotted 
and the solid curves in the bottom panels. One also ob- 
serves that the distribution depends on whether the mass 
(symmetry energy) evolves during the evaporation or not. 
The SMM calculation with the standard de-excitation (i.e 
the old de-excitation) leads to a narrow final distribution 
and the isotopes are concentrated close to the /^-stability 
line. The difference in the final yield distributions for 7 = 
15 MeV and 7 = 25 MeV is very small. This difference is 
however much more pronounced in the new de-excitation 
calculation. The final isotope distributions in this case 
are considerably wider, and shifted toward neutron-rich 
side. The SMM calculation assuming the mass (symme- 
try energy) evolution during the evaporation therefore 
leads to larger yields for neutron rich fragments. A simi- 
lar observation was made by Buyukcizmeci et al. j63| in 
their calculation of the primary and secondary fragment 
isotopic distributions in 197 Au, 124 Sn and 124 La systems. 
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By using the experimental masses during the evaporation 
from the primary fragments, the emission of the charged 
particles are suppressed by the binding energy and the 
Coulomb barrier. In the case of small 7, the binding en- 
ergy in the beginning of evaporation process essentially 
favors emission of charged particles. When the nucleus 
has cooled down sufficiently to restore the normal sym- 
metry energy, the remaining excitation energy is rather 
low (E* I A < 1 MeV) for the nucleus to evaporate many 
neutrons. 

The above comparison of the experimentally observed 
isoscaling properties with the statistical multifragmen- 
tation model shows that, irrespective of the secondary 
de-excitation, the final fragment distribution depends 
strongly on the available free energy and the strength 
of the symmetry energy. A significantly lower value of 
the symmetry energy than that assumed for cold isolated 
nuclei is required to explain the isotopic compositions of 
the fragments produced in multifragmcntation reaction. 
The difference between the two kinds of evaporation cal- 
culations gives a measure of the uncertainty expected in 
the present analysis. The results above indicate that the 
properties of nuclei produced at high excitation energy, 
isospin and reduced density could be significantly differ- 
ent from those of the cold isolated nuclei. Such informa- 
tion can provide important inputs for the understanding 
of the nuclear composition of supernova matter where hot 
and neutron rich nuclei are routinely produced 0, Il2j . 

VI. SUMMARY AND CONCLUSIONS 

In summary, we have measured the isotopic yield dis- 
tribution of the fragments produced in the multifragmen- 



[1] D.V. Shetty, A.S. Botvina, S.J. Yennello, A. Keksis, E. 
Martin and G.A. Souliotis, Phys. Rev. C 71, 024602 
(2005). 

[2] A. LeFevre et at, Phys. Rev. Lett. 94, 162701 (2005). 
[3] D. Henzlova, A.S. Botvina, K.H. Schmidt, V. Henzl, P. 

Napolitani, M.V. Ricciardi, nucl-ex/0507003 (2005). 
[4] A.S. Botvina and I.N. Mishustin, Phys. Lett. B 584, 233 

(2004). 

[5] A.S. Botvina, O.V. Lozhkin, and W. Trautmann, Phys. 

Rev. C 65, 044610 (2002). 
[6] J. Lattimer, C. Pethick, M. Prakash and P. Hansel, Phys. 

Rev. Lett. 66, 2701 (1991). 
[7] C. Lee, Phys. Rep. 275, 255 (1996). 
[8] C.J. Pethick and D.G. Ravenhall, Annu. Rev. Nucl. Part. 

Sci. 45, 429 (1995). 
[9] J.M. Lattimer and M. Prakash, Phys. Rep. 333, 121 

(2000) . 

[10] W.R. Hix et at, Phys. Rev. Lett. 91, 201102 (2003). 
[11] J.M. Lattimer and M. Prakash, Astrophys. J. 550, 426 

(2001) . 

[12] A.S. Botvina and I.N. Mishustin, Phys. Rev. C 72, 
048801 (2005). 



tation of 40 Ar, 40 Ca + 58 Fe, 58 Ni reactions at 25 - 53 
MeV/nuclcon. The symmetry energy and the isoscal- 
ing properties of the fragments produced were studied 
within the framework of statistical multifragmentation 
model. It is observed that the isoscaling parameter a 
for the hot fragments decrease with increasing excitation 
energy and decreasing symmetry energy. The a values 
increase with increasing difference in the isospin of the 
fragmenting system. Similar behavior is also observed 
for the cold secondary fragments. The sequential decay 
of the primary fragments to secondary fragments is ob- 
served to have very little influence on the isoscaling pa- 
rameter as a function of excitation energies and isospin 
of the fragmenting system. The symmetry energy how- 
ever, strongly influences the isospin properties of the hot 
fragments. The experimentally determined scaling pa- 
rameters could be explained by a symmetry energy that 
is as low as 10 - 15 MeV, and significantly lower than 
that for the normal (cold) nuclei at saturation density. 
The present results indicate that the isospin properties 
of the fragments produced at high excitation energy and 
reduced density in multifragmentation reaction are sen- 
sitive to the symmetry energy. 

VII. ACKNOWLEDGMENT 



The authors wish to thank the Cyclotron Institute 
staff for the excellent beam quality. This work was 
supported in parts by the Robert A. Welch Foundation 
through grant No. A-1266, and the Department of En- 
ergy through grant No. DE-FG03-93ER40773. 



[13] B.A. Brown, Phys. Rev. Lett 85, 5296 (2000). 
[14] C.J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 
5647 (2001). 

[15] R.J. Furnstahl, Nucl. Phys. A 706, 85 (2002). 

[16] K. Oyamatsu, I. Tanihata, Y. Sugahara, K. Sumiyoshi 

and H. Toki, Nucl. Phys. A 634, 3 (1998). 
[17] A. Ono, P. Danielewicz, W.A. Friedman, W.G. Lynch 

and M.B. Tsang, Phys. Rev. C 68, 051601 (2003). 
[18] L.W. Chen, CM. Ko and B.A. Li, Phys. Rev. Lett. 94, 

032701 (2005). 
[19] D.V. Shetty et at, Phys. Rev. C 70, 011601 (2004). 
[20] M.B. Tsang et ai, Phys. Rev. Lett. 92, 062701 (2004). 
[21] B.A. Li, CM. Ko and W. Bauer, Int. J. Mod. Phys. E 7, 

147 (1998). 

[22] V. Baran, M. Colonna, V. Greco and M. DiToro, Phys. 

Rep. 410, 335 (2005). 
[23] V. Baran, M. Colonna, M. Di Toro, V. Greco, M. 

Zielinska-Pfabe and H.H. Wolter, Nucl. Phys. A 703, 603 

(2002). 

[24] T.X. Liu et ai, Phys. Rev. C 69, 014603 (2004). 
[25] M.B. Tsang et ai, Phys. Rev. C 64, 054615 (2001). 
[26] M.B. Tsang et ai, Phys. Rev. Lett. 86, 5023 (2001). 



14 



[27] W.P. Tan et al, Phys. Rev. C 64, 051901 (2001). 

[28] D. Lunney, J.M. Pearson and C. Thibault, Rev. Mod. 

Phys. 75, 1021 (2003). 
[29] A.E.L. Dieperink, Y. Dewulf, D. Van Neck, M. Waroquicr 

and V. Rodin, Phys. Rev. C 68, 064307 (2003). 
[30] W. Zuo, I. Bombaci and U. Lombardo, Phys. Rev. C 60, 

024605 (1999). 

[31] M. Brack, C. Guet and H.B. Hakansson, Phys. Rep. 123, 
276 (1985). 

[32] J.M. Pearson and R.C. Nayak, Nucl. Phys. A 668, 163 
(2000) 

[33] W.D. Myers and W.J. Swiatecki, Nucl. Phys. 81, 1 
(1966). 

[34] K. Pomorski and J. Dudek, Phys. Rev. C 67, 044316 
(2003). 

[35] B. Borderie et al, Phys. Rev. Lett. 86, 3252 (2001). 
[36] D.R. Bowman et al, Phys. Rev. Lett. 67, 1527 (1991). 
[37] M. D'Agostino et al, Phys. Lett. B 371, 175 (1996). 
[38] L. Beaulieu et al, Phys. Rev. Lett. 84, 5971 (2000). 
[39] W.P. Tan, S.R. Souza, R.J. Charity, R. Donangelo, 

W.G. Lynch, and M.B. Tsang, Phys. Rev. C 68, 034609 

(2003). 

[40] H. Xi, W.G. Lynch, M.B. Tsang, W.A. Friedman, and 
D. Durand, Phys. Rev. C 59, 1567 (1999). 

[41] H. Xi, W.G. Lynch, M.B. Tsang, and W.A. Friedman, 
Phys. Rev. C 54, 2163 (1996). 

[42] H. Johnston et al, Phys. Rev. C 56, 1972 (1997). 

[43] R.P. Schmitt et al, Nucl. Instrum. Methods Phys. Res. 
A 354, 487 (1995). 

[44] Ad.R. Raduta, Eur. Phys. Jour. A24, 85 (2005). 



[45] M.B. Tsang et al, Phys. Rev. C 64, 041603 (2001). 

[46] S. Albergo et al, Nuovo Cimento A89, 1 (1985). 

[47] J.P. Bondorf, A.S. Botvina, A.S. Iljinov, I.N. Mishustin 

and K. Sneppen, Phys. Rep. 257, 133 (1995). 
[48] J. Randrup and S. Koonin, Nucl. Phys. A 356, 223 

(1981). 

[49] A.S. Botvina, A.S. Iljinov, I.N. Mishustin, J.P. Bondorf, 
R. Donangelo and K. Sneppen, Nucl. Phys. A 475, 663 
(1987). 

[50] D.H.E. Gross et al, Rep. Progr. Phys. 53, 605 (1990). 

[51] D.H.E. Gross et al, Z. Phys. A 309, 41 (1982). 

[52] X.Z. Zhang, D.H.E. Gross, S. Xu, and Y.M. Zheng, Nucl. 
Phys. A 461, 641 (1987). 

[53] Isospin Physics in Heavy Ion Collisions at Intermediate 
Energies, edited by B.-A. Li and W. Schroder (Nova Sci- 
ence, New York, 2001). 

[54] P. Chomaz and F. Gulminelli, Phys. Lett. B 447, 221 
(1999). 

[55] A.S. Botvina et al, Nucl. Phys. A 584, 737 (1995). 
[56] R.P. Scharenberg et al, Phys. Rev. C 64, 054602 (2001). 
[57] N. Bellaize et al, Nucl. Phys. A 709, 367 (2002). 
[58] S.P. Avdeyev et al, Nucl. Phys. A 709, 392 (2002). 
[59] J. Pochodzalla et al, Phys. Rev. Lett. 75, 1040 (1995). 
[60] M. D'Agostino et al, Nucl. Phys. A 650, 329 (1999). 
[61] D. Cussol et al, Nucl. Phys. A 561, 298 (1993). 
[62] A.S. Botvina and I.N. Mishustin, Phys. Rev. C 63, 
061601 (2001). 

[63] N. Buyukcizmeci, R. Ogul, and A.S. Botvina, Eur. Phys. 
J. 25, 57 (2005). 



